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Inflation 

1.1. Tkachev 0, 

a CERN Theory Division, CH-1211 Geneva 23, Switzerland. 

Inflationary cosmology is reviewed. Particular attention is given to processes of creation. 


1. Introduction 

A typical universe should have had Plankian 
size, live Plankian time and contain 1 parti¬ 
cle. Yet, the observable Universe contains 10 90 
particles in it and had survived 10 65 Plankian 
times. Where does it all came from ? In other 
words, why is the Universe so big, flat and old 
(t > 10 10 years), homogeneous and isotropic 
(ST/T ~ 10 -5 ), why does it contain so much 
entropy and does not contain unwanted relics? 
These puzzles of classical cosmology were solved 
with invention of Inflation [Q. 

1.1. Getting something for nothing 

Stress-energy tensor T which drives the ex¬ 
pansion of a homogeneous universe can be charac¬ 
terised by two parameters, energy density, T 9 = 
p, and pressure T'- = — pS *•. Conservation of 
energy and momentum T^ v , v = 0 in an ex¬ 
panding Friedmann universe takes a simple form 
p + 3 H(p + p) = 0. Here H is a Hubble pa¬ 
rameter, H = a/a , and a(t) is a scale factor 
which describes the expansion of the Universe. 
Consider the stress-energy tensor T of a vac¬ 
uum. Vacuum has to be Lorentz invariant, hence 
T" = V 6" and we find p = —p. Therefore, the 
energy of a vacuum stays constant despite the 
expansion. In this way, room for matter full of 
energy could have being created. It remains to 
find out how vacuum energy was converted into 
radiation at some later stage. 

1.2. Horizon problem and the solution 

The size of a causally connected region (hori¬ 
zon) scales in proportion to time, 1 ?h oc t. On 
the other hand, the physical size of a given patch 
grows in proportion to the scale factor, Rp oc 
a{t) oc f 7 . Exponent 7 depends upon equation 


of state, 7 = 1/2 for radiation and 7 = 2/3 for 
matter dominated expansion. In any case, for the 
“classical” Friedmann Universe 7 < 1 and hori¬ 
zon expands faster than volume. Take the largest 
visible patch today. It follows that in the past 
it should have contained many casually discon¬ 
nected regions and the question arises why the 
Universe is so homogeneous at large scales ? This 
problem can be solved if during some period of 
time the volume had expanded faster than the 
horizon. During such a period the whole visi¬ 
ble Universe can be inflated from one (“small”) 
causally connected region. Clearly, this happens 
if 7 > 1, which means a > 0. Either of these two 
conditions can be used as definition of an infla¬ 
tionary regime. 

Using the Einstein equation d = —47rGa(p + 
3p)/3 we find that the inflationary stage is real¬ 
ized when p < —p/3. If p = —p we have de Sitter 
metric and the Universe expands exponentially. 

A crucial and testable prediction of Inflationary 
cosmology is a flat Universe, H = 1 (as usual, fl(f) 
is the ratio of current and critical energy densities, 
n(i) = 8nGp/3H 2 ). Indeed, Einstein equations 
can be cast into the form d 2 (fl(f) — 1) = do(flo — 
1). Accelerated expansion, a > 0, increases a and 
therefore drives fl(t) to 1. 

The first model of inflation was de-facto sug¬ 
gested in ||. De Sitter expansion appeared as a 
result of vacuum polarization effects in a one-loop 
order - too complicated to be sure that higher or¬ 
der corrections are unimportant. 

1.3. Arranging for a vacuum 

Consider X/„ of a real scalar field cf> 

T 
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with the Lagrangian C = ^d^d^cj) — V (<{>). In 
a state when all derivatives of (j) are negligible, 
d^cj) ~ 0 , the stress-energy tensor of a scalar field 
is that of a vacuum, ~ V(^>) < 7 ^ . 

There are two basic ways to arrange 4> — const 
and hence to imitate the vacuum-like state. 

1. The simplest possibility was suggested by 
A. Guth in his original paper |l|. Consider po¬ 
tential V (</>) which has a local minimum with a 
non-zero energy density separated from the true 
ground state by a potential barrier. The Uni¬ 
verse will be trapped in the meta-stable minimum 
for a while and expansion will diminish all field 
gradients. Then the Universe enters a vacuum 
state. Subsequent phase transition into the true 
minimum ends inflation and creates the radiation 
phase. Today the model of Guth and its variants 
based on potential barriers is good for illustration 
purposes only. It did not stand up to observations 
since inhomogeneities which are created during 
the phase transition into the radiation phase are 
too large ||. 

2. A. Linde was first to realise that things work 
in the simplest possible set up |Q. Consider po¬ 
tential 

= \mltf . ( 1 ) 

Equation of the field motion in expanding Uni¬ 
verse is (j> + 3 H(j) + m^</> = 0 . If H » to the 
“friction” is too big and the field (almost) does 
not move. Therefore time derivatives in can 
be neglected and inflation starts (in sufficiently 
homogeneous patch of the Universe). Hubble pa¬ 
rameter in this case is given by H « m<f>/Mp\ and 
we see that inflation starts if initial field value 
happen to satisfy </> > Mp\. During inflationary 
stage the field slowly rolls down the potential hill. 
This motion is very important in the theory of 
structure creation. Inflation ends when <f> ~ Mp\. 
At this time field oscillations start around poten¬ 
tial minimum and latter decay into radiation. In 
this way matter was likely created in our Uni¬ 
verse. 

2. Unified theory of creation 

Small fluctuations of any field obey 
Uk + [fc 2 + Wgfj(T)] Uk = 0 . ( 2 ) 


Effective mass is time dependent here because of 
the expansion of the Universe. Because m e ff is 
time dependent, it is not possible to keep fluctu¬ 
ations in vacuum. If one arranges to put oscil¬ 
lators with momentum k into the vacuum, they 
will not be in vacuum at a latter time since this 
vacuum would correspond to wrong value of the 
field mass. 

Some remarks are in order: 

• Eq. (|j) is valid for all particle species 

• Equation looks that simple in conformal 
reference frame ds 2 = a(r) 2 ( dr 2 — dx 2 ). 
Everywhere in this chapter a “dot” means 
derivative with respect to r. 

• Of particular interest are ripples of space- 
time itself 

— curvature fluctuations (scalar fluctua¬ 
tions of the metric) 

— gravitons (tensor fluctuations of the 
metric) 

• TO e ff may be non-zero even for massless 
fields 

— graviton is the simplest example ||, 

m eff = ~ d / a 

— m 2 ff for curvature fluctuations has 
similar structure || with a being re¬ 
placed by acj)/H 

• For conformally coupled, but massive scalar 
TOeff = Too a(r) 

Creation was only possible because nature is not 
conformally-invariant. Otherwise m e s = 0 and 
vacuum remains vacuum forever. 

2.1. Sources of creation 

Amplitudes Uk in Eq. (J2j) are quantum opera¬ 
tors and a theory of creation reduces to the the¬ 
ory of Bogolyubov transformations or to a theory 
of particle creation in homogeneous time varying 
classical background. There are two important 
instances of such background in cosmology: 

• Expansion of space-time, a(r) 
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• Motion of the inflaton field, </>(t) 

Both can be operational at any epoch of creation 

• During inflation 

• While the inflaton oscillates (reheating) 

During inflation superhorizon size perturbations 
of metric are created which give seeds for Large 
Scale Structure (LSS) formation and eventually 
lead to formation of galaxies, the Solar system 
and all the rest which we can see around us. 
During reheating matter itself is created. Over¬ 
all there are four different situations (two sources 
times two epochs). If coupling to the inflaton is 
not essential, the corresponding process will be 
called “pure gravitational creation” in what fol¬ 
lows. Let me consider all four possibilities in turn, 
starting from 


P(k) T = 


16 H 2 

TV Mpj 


( 4 ) 


These spectra can be approximated as power law 
functions 


P(k)s = P(k 0 ) s 


P(k)r = P(k 0 )T 



( 5 ) 


where scale dependence k enters via weak depen¬ 
dence of H on the current field value (current 
means here corresponding to the moment when 
the scale gets bigger than the horizon and evo¬ 
lution of the mode k freezes out). One expects 
spectra to be nearly scale invariant since the field 
(almost) does not move. Indeed, expanding in 
slow roll parameters one finds 


2.2. Gravitational creation of metric per¬ 
turbations [7j] 

During inflation the motion of the inflaton field 
is slow, while the expansion of the Universe is 
fast. It follows that relevant cosmological scales 
encompass small A <f> interval. E.g. in the model 
Eq. ( 0 ) the whole visible Universe is inflated away 
while the inflaton field <p changes from 4Mpi to 
Mpi. Potential V(cf>) should be relatively “flat” 
over this range of A<p to maintain the inflationary 
regime. We may conclude that observables should 
essentially depend on a first few derivatives of V 
and the shape of the potential outside this region 
of A (f> is irrelevant. One may construct dimen¬ 
sionless quantities (slow roll parameters) out of 
potential derivatives 


e _ M 2 X 

(V 
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V" 

V ~ 87T 

V ' 


The value of potential itself during this period is 
also relevant and defines the value of the Hubble 
parameter, H 2 (4>) = 87 rGU((^)/ 3 . 

Solutions of Eq. (|2|) with vacuum initial condi¬ 
tions give for the power spectra of scalar (curva¬ 
ture) and tensor (gravity waves) perturbations 


P(k) s 


1 H 2 
ne Mp j ’ 


n — 1 = 2rj — 6e , 

ut = — 2e . (6) 

2.2.1. Consistency relation 

According to Eqs. ([|) the ratio of power in 
tensor to scalar perturbations is equal to 16e. 
On the other hand the exponent of tensor per¬ 
turbations is also proportional to e, see Eq. ( 0 ). 
This gives “consistency relation”, tit — — 0.14r, 
where the ratio of tensor to scalar power is ex¬ 
pressed through the ratio of directly measurable 
respective contributions to quadrupole CMBR 
anisotropy, r = C T /C s . Verification of the con¬ 
sistency relation should give major peace of evi¬ 
dence that the Inflation did happened. 

2.2.2. Testing Inflation 

Predictions of inflationary theory can be tested 
measuring CMBR anisotropy and power spectra 
of galaxies distribution. All tests completed so 
far are in agreement with predictions. Latest 
CMBR data [0 imply that the Universe is flat 
Do = 1.02 ± 0.04. The spectrum of perturbations 
is nearly scale invariant, n s = 0.97 ±0.1, and 
Gaussian |)j. Overall normalisation of CMBR 
spectrum fixes inflaton parameters. In the sim¬ 
plest case of Eq. (0) m$ ~ 10 13 GeV. 

Consistency relation was not tested yet since 
it requires measurement of tensor perturbations. 
In forthcoming CMBR experiments tensor and 
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scalar modes may be disentangled. This can be 
done unambiguously because tensor and scalar 
perturbations contribute differently to the polar¬ 
isation of CMBR (for a review see @) which 
hopefully will be measured. 

Measured position of the first acoustic peak in 
CMBR not only tells us that the Universe is flat, 
but also that the isocurvature perturbations are 
ruled out as the primary block of structure forma¬ 
tion. While during inflation both types of pertur¬ 
bations, curvature and isocurvature, can be pro¬ 
duced, the models which do not involve inflation 
(e.g. network of cosmic strings) produce isocur¬ 
vature fluctuations. Therefore this is an impor¬ 
tant test of inflation ]fl| . In principle, curvature 
perturbations may be mimicked by causal pro¬ 
cesses 0 , but polarisation measurements will 
also give unambiguous proof that density pertur¬ 
bations are of superhorazon origin [ p~3[ ■ 

While the tensor mode is not measured yet, it 
is restricted by CMBR anisotropy measurements. 
Different models of inflation occupy well defined 
and different regions in the (r,n) parameter plane 

[14] , see Eqs. (j|), (Q). Regions, specified for new 

[15] , chaotic ||] and hybrid [01 inflationary mod¬ 
els are shown in Fig. [I], Other inflationary models 
are also bounded to these regions. For example, 
“natural” inflation 0 ] fits the region of “new” 
inflationary model. Parameter space favoured by 
current CMBR and LSS measurements is also 
shown [0 as a shaded area (note however that 
size and shape of this region depends upon pri¬ 
ors used). E.g., for the V oc (jP chaotic inflation 
model all this means that p < 6 at 99% c.l. |l 8 | , 
while the hybrid model of inflation is disfavoured. 

2.3. Gravitational creation of matter 

Let us consider the gravitational creation of 
matter. The source of creation is m e g = mo a(r). 
While matter is created in tiny amounts by this 
process, it is most effective for heaviest particles. 
If such superheavy particles do exist in nature, 
this process naturally leads to Superheavy Dark 
Matter (SDM) 0,0. If these particles are un¬ 
stable but long lived, they can explain [0 puz¬ 
zling Ultra-High Energy Cosmic Ray (UHECR) 
events. 



0.6 0.8 1.0 1.2 


TX 

Figure 1. Different models of inflation pre¬ 
dict n and r to be found in specified regions. 
Shaded area indicates parameter range favoured 
by CMBR and LSS measurements. 


2.3.1. Friedmann cosmology 

It is particle mass which couples quantum fluc¬ 
tuations to the background expansion and serves 
as the source of particle creation. Therefore we 
expect for the abundance of Super-Heavy parti¬ 
cles 

n S H oc m| H a -3 . (7) 

It is the expansion of the Universe which causes 
particle creation, therefore creation is most effi¬ 
cient at time tq when H « tosh, while (comov¬ 
ing) particle number is an adiabatic invariant at 
later times. We expect that coefficient in Eq. 
(j7j) should not be much smaller unity if we nor¬ 
malize a(ro) = 1. It follows that stable particles 
with tosh > 10 9 GeV will overdose the Universe 
if there was Friedmann singularity in the past - 
Friedmann cosmology and SDM mutually exclude 
each other ( 0 . 

2.3.2. Inflationary cosmology 

In inflationary cosmology there is no singu¬ 
larity and the Hubble constant is limited, H <, 
m Therefore, the production of particles with 
msH > H ~ 10 13 GeV is suppressed. 

The present day ratio of the energy density 
in S'.ff-particles to the critical energy density is 
shown in Fig. |] [0 (see also 0). Super-Heavy 
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Figure 2. Present day density of gravitation¬ 
ally created Superheavy Dark Matter. Solid 
lines correspond to minimally or conformally cou¬ 
pled bosons, dotted line describes abundance of 
fermions. 


particles with the mass few x 10 13 GeV are ex¬ 
cellent candidates for SDM and progenitors of 
UHECR. If it will be proven that UHECR are due 
to decays of SDM, it will mean that Friedmann 
expansion was preceded by some other epoch, 
likely by Inflation. 

Consistency with observations requires 
flsuh 2 ^ 0.3. We see that light particles, 
tosh ■C 10 13 GeV, are overproduced by many 
orders of magnitude, unless their coupling to 
curvature is conformal £ = 1/6. This poses 
a serious danger since supresymmetry and 
suprgravity models predict many such particles, 
e.g. moduli, gravitino, etc. They may be dan¬ 
gerous relics, even if unstable, since they can 
easily survive to post-nucleosinthesis epoch and 
their decay products destroy 4 He and D nuclei 
by photodissociation 1231. 


2.4. Decay of the inflaton oscillations 

While bosons and fermions are created equally 
efficiently by a pure gravitational mechanism, 
coupling to inflaton uncovers deep differences be¬ 
tween them. Effective mass of a scalar A' (in¬ 
teraction Lagrangian L lnt = | g 2 (j> 2 X 2 ) and of a 
fermion (interaction Lagrangian L; Ilt = 


is given by the following expressions 

scalarA' : m 2 ff = m 2 x + g 2 (j> 2 {t) (8) 

fermion : m e s = m^ + g<p(t ). (9) 


Effective mass of a scalar A depends quadrati- 
cally upon the inflaton field strength and there¬ 
fore it is always larger than the bare mass mx ■ In 
the case of fermions inflaton field strength enters 
linearly and effective mass can cross zero. Even 
superheavy fermions can be created easily during 
these moments of zero crossing [Q. (It is eas¬ 
ier to create a light field and it is effective mass 
which counts at creation. At the end of the day 
it is bare mass which counts.) 

Coupling g by itself is not relevant for the pro¬ 
cess of creation, g always comes in combination 
with inflaton field strength. To make dimension¬ 
less combination out of it we have to re-scale gcf> 
by a typical time scale relevant for creation. In 
the present case this will be period of inflaton os¬ 
cillations or inverse inflaton mass 
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g = 


g 2 4> 2 

4 ml 


( 10 ) 


Parameter q determines the strength of particle 
production caused by the oscillations of the infla¬ 
ton field. It can be very large ||25|| even when g is 
small since </> 2 /m 2 ss 10 12 . 


2.4.1. Matter creation: Bose versus Fermi 

Bose-stimulation aids the process of creation 
of bosons. Occupation numbers grow exponen¬ 
tially with time, n = e Mt , which results in a 
fast, explosive decay of inflaton 25 2tJ and cre¬ 
ates large classical fluctuations of all Bose-fields 
involved. This can have a number of observ¬ 
able consequences: non-thermal phase transitions 
J27],|2^], generation of a stochastic background of 
the gravitational waves 291, and a possibility for 
a novel mechanism of baryogenesis 30 are some 
examples. 

Fortunately, the system in this regime of parti¬ 
cle creation became classical and can be studied 
on a lattice |il] HH . 

On the other hand in the case of fermions n < 1 
at all times because of the Pauli blocking. This 
may create an impression that the fermionic chan¬ 
nel of inflaton decay is not important. This is a 
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Figure 3. Effectiveness of Super-Heavy particle 
production. Solid lines production of fermions, 
dotted lines production of bosons. 


false impression |34| . Production of fermions can 
be more efficient compared to bosons J 24 J. Van¬ 
ishing effective mass allows for a larger Fermi- 
momentum of created particles and fermions can 
outnumber bosons. 

A fraction of the initial energy density which 
goes to bosons and fermions respectively is plot¬ 
ted in Fig. | as a function of particle mass for 
several values of q. Indeed, superheavy fermions 
are more efficiently created compared to bosons 
at the same value of q. 

2.4.2. Non-thermal phase transitions 

Bosons of moderate or small mass (compared to 
inflaton), and sufficiently large coupling (q > 10 4 
or g 2 > 10 -8 Q) are produced explosively, 
their density grows exponentially. Field vari¬ 
ances, ( X 2 ), can reach large values, larger than 
in a thermal equilibrium at the same energy den¬ 
sity. These were calculated for different values of 
parameters in Ref. fu| . For one choice of mx 
and q the time dependence of field variances is 
shown in Fig. [|. For smaller mx and q variances 
can reach larger values. By the final moment of 
time shown in this figure inflaton zero mode had 
already decayed, but the system is still very far 
from thermal equilibrium. 

In the theory of phase transitions the effective 
mass of the Higgs field in a medium is = 
—ft 2 + h 2 (X 2 ). If ( X 2 ) is sufficiently large a phase 



Figure 4. Late-time dependence of field variances 
and of inflaton zero mode </>o- Time and mx are 
in units of m^. At r = 0 variances are in vacuum. 


transition can occur. During preheating in the 
process of inflaton decay this can have especially 
interesting consequences. Symmetry restoration 
and subsequent breaking with possible formation 
of topological defects j|(| can happen even in 
GUTs with large symmetry breaking scale. 

2.5. Creation of matter during Inflation 

While inflaton during Inflation moves slowly, 
it does move. This motion may lead to the cre¬ 
ation of particles coupled to it. Most effectively 
and with observable consequences this occurs in 
the case of superheavy fermions. Features in the 
power spectrum of perturbations which may leave 
trace in CMBR and LSS appear at scales corre¬ 
sponding to zero crossings of effective mass. This 
can be a probe of Sub-Plankian particle content. 
Multiplet of N particles with ~ Mpi and 
coupling g > 0.2/N 2 / 5 is detectable J|7|. 

3. Conclusions 

Inflationary cosmology is a beatuful theory. It 
inputs unknown and (largely) arbitrary initial 
conditions and replaces them with testable pre¬ 
dictions. Even the simplest inflationary model, 
Eq. (lj), reproduces all relevant features of ob¬ 
servable Universe after adjustment of a single pa¬ 
rameter, ?Ti 0 ~ 10 13 GeV. 









7 


REFERENCES 


1 . 

2 . 

3. 

4. 

5. 

6 . 
7. 


A. H. Guth, Phys. Rev. D 23 (1981) 347. 

A. A. Starobinsky, JETP Lett. 30 (1979) 682; 
Phys. Lett. B 91 (1980) 99. 

A. H. Guth and E. J. Weinberg, Phys. Rev. 
D 23 (1981) 876. 

A. D. Linde, Phys. Lett. B 129 (1983) 177. 
L. P. Grishchuk, Sov. Phys. JETP 40 (1975) 
409. 

V.N. Lukash, Sov. Phys. JETP 52 (1980) 807. 
For a reviews and numerous references see: 
V. F. Mukhanov, H. A. Feldman and 
R. H. Brandenberger, Phys. Rept. 215 (1992) 
203; J. E. Lidsey et a/., Rev. Mod. Phys. 69 
(1997) 373; D. H. Lyth and A. Riotto, Phys. 


8 . C. B. Netterfield et al. [astro- 

ph/010446C; 

C. Prykeet al. astro-ph/0104490 

; S. Hanany 


et al ., Ap. J. 545 (2000) L5; A. H. Jaffe et al. 
Phys. Rev. Lett. 86 (2001) 3475. 

9. J. H. Wu et al, |astro-ph/0104248| . 

10. W. Hu and M. J. White, New Astron. 2 
(1997) 323. 

11. R. G. Crittenden and N. Turok, Phys. Rev. 
Lett. 75 (1995) 2642; A. Albrecht et al. 
Phys. Rev. Lett. 76 (1996) 1413; W. Hu and 
M. J. White, Phys. Rev. Lett. 77 (1996) 1687. 

12. N. Turok, Phys. Rev. Lett. 77 (1996) 4138. 

13. D. N. Spergel and M. Zaldarriaga, Phys. Rev. 
Lett. 79 (1997) 2180. 

14. S. Dodelson, W. H. Kinney and E. W. Kolb, 
Phys. Rev. D 56 (1997) 3207; see also 
V. N. Lukash and E. V. Mikheeva, Int. J. 
Mod. Phys. A 15 (2000) 3783. 

15. A. D. Linde, Phys. Lett. B 108 (1982) 389; 
A. Albrecht and P. J. Steinhardt, Phys. Rev. 
Lett. 48 (1982) 1220. 

16. A. D. Linde, Phys. Rev. D 49 (1994) 748. 

17. K. Freese, J. A. Frieman and A. V. Olinto, 
Phys. Rev. Lett. 65 (1990) 3233. 

18. W. H. Kinney, A. Melchiorri and A. Riotto, 
Phys Rev. D 63 (2001) 023505; S. Hannestad 


et al. astro-ph/0103047 . 


19. D. J. Chung, E. W. Kolb and A. Riotto, Phys. 
Rev. D 59 (1999) 023501. 

20. V. Kuzmin and I. Tkachev, JETP Lett. 68 
(1998) 271; Phys. Rev. D 59 (1999) 123006. 


21. V. Berezinsky, M. Kachelriess and 
A. Vilenkin, Phys. Rev. Lett. 79 (1997) 
4302; V. A. Kuzmin and V. A. Rubakov, 
Phys. Atom. Nucl. 61 (1998) 1028. 

22. R. Kallosh, L. Kofrnan, A. D. Linde and 
A. Van Proeyen, Phys. Rev. D 61 (2000) 
103503; G. F. Giudice, I. Tkachev and A. Ri¬ 
otto, JHEP 9908 (1999) 009; H. P. Nilles, 
M. Peloso and L. Sorbo, Phys. Rev. Lett. 87 
(2001) 051302 

23. J. R. Ellis, A. D. Linde and D. V. Nanopou- 
los, Phys. Lett. B 118 (1982) 59; D. Lind- 
ley, Astrophys. J. 294 (1985) 1; J. R. Ellis, 
D. V. Nanopoulos and S. Sarkar, Nucl. Phys. 
B 259 (1985) 175. 

24. G. F. Giudice et al. JHEP 9908 (1999) 014. 

25. L. Kofrnan, A. D. Linde and A. A. Starobin¬ 
sky, Phys. Rev. Lett. 73 (1994) 3195; Phys. 
Rev. D 56 (1997) 3258. 

26. Y. Shtanov, J. Traschen and R. H. Branden¬ 
berger, Phys. Rev. D 51 (1995) 5438. 

27. L. Kofrnan, A. D. Linde and A. A. Starobin¬ 
sky, Phys. Rev. Lett. 76 (1996) 1011. 

28. I. I. Tkachev, Phys. Lett. B 376 (1996) 35. 

29. S. Y. Khlebnikov and I. I. Tkachev, Phys. 
Rev. D 56 (1997) 653. 

30. E. W. Kolb, A. D. Linde and A. Riotto, Phys. 
Rev. Lett. 77 (1996) 4290; E. W. Kolb, A. Ri¬ 
otto and I. I. Tkachev, Phys. Lett. B 423 
(1998) 348; J. Garcia-Bellido et al, Phys. 
Rev. D 60 (1999) 12 3504; J. Garcia-B ellido 
and E. Ruiz Morales, hep-ph/0109230| . 

31. S. Y. Khlebnikov and I. I. Tkachev, Phys. 
Rev. Lett. 77 (1996) 219; Phys. Rev. Lett. 
79 (1997) 1607. 

32. T. Prokopec and T. G. Roos, Phys. Rev. D 
55 (1997) 3768. 

33. G. N. Felder and I. Tkachev, hep-ph/0011159| . 

34. P. B. Greene and L. Kofrnan, Phys. Lett. B 
448 (1999) 6 ; J. Baacke, K. Heitmann and 
C. Patzold, Phys. Rev. D 58, 125013 (1998). 

35. S. Y. Khlebnikov and I. I. Tkachev, Phys. 
Lett. B 390 (1997) 80. 

36. I. Tkachev et al, Phys. Lett. B 440 (1998) 
262; S. Khlebnikov et al, Phys. Rev. Lett. 81 
(1998) 2012. 

37. D. J. Chung et al, Phys. Rev. D 62 (2000) 
043508. 

















